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Abstract 

Structural optimization methods in MSC /NASTRAN 
are used to size substructures and to reduce the weight 
of a composite sandwich cryogenic tank for future 
launch vehicles. Because the feasible design space of 
this problem is non-con vex, many local minima are 
found. This non-convex problem is investigated in 
detail by conducting a series of analyses along a design 
line connecting two feasible designs. Strain constraint 
violations occur for some design points along the 
design line. Since MSC/NASTRAN uses gradient- 
based optimization procedures, it does not guarantee 
that the lowest weight design can be found. In this 
study, a simple procedure is introduced to create a new 
starling point based on design variable values from 
previous optimization analyses. Optimization analysis 
using this new starting point can produce a lower 
weight design. Detailed inputs for setting up the 
MSC/NASTRAN optimization analysis and final tank 
design results are presented in this paper. Approaches 
for obtaining further weight reductions are also 
discussed. 

Introduction 

The Held of structural optimization is well developed 
and is increasingly used in automotive, aerospace, 
machine design, civil, and other engineering fields 
1 1,2,3). Structural optimization methods are useful for 
conducting trade-off studies and sizing of structures. 
Future launch vehicles must be highly reliable 
lightweight structures to achieve the goal of affordable 
space travel. Designers of future launch vehicles need 
to take full advantage of optimization methods to obtain 
lightweight vehicle designs and evaluate them for 
meeting high reliability design requirements. 

In future launch vehicles, cryogenic tanks will not only 
carry liquid oxygen (LOX) and hydrogen (LH2) fuel, 
but also function as the primary vehicle load carrying 
members. Current designs of lifting body type launch 
vehicles use multi-lobed tank configurations [4J. Such 
designs require the use of external rib-like standoff 
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structures to support thermal protection tiles over the 
lobed intersections to create smooth aerodynamic 
surfaces. A recent study at NASA Langley Research 
Center (LaRC) found that the fuel volume efficiency 
could he enhanced by developing a conformal tank that 
closely follows the outside mold line of the vehicle [4 1 . 
The LaRC's conformal cryogenic tank design and its 
dimensions are shown in Figure 1. 

The objective of this study is to use the optimization 
capabilities available in MSC/NASTRAN [5| for sizing 
and weight reduction of honeycomb sandwich 
conformal cryogenic tank structures. The facesheet and 
honeycomb core thicknesses of the sandwich structures 
are used as design variables to be optimized. The 
NAvSTRAN input entries for optimizing these design 
variables are developed and tested. Technical 
difficulties encountered in using optimization for sizing 
and finding optimal low weight design of the conformal 
tank arc discussed. Nonconvexity of the feasible design 
space, which can often have local minima, is 
investigated. A simple procedure for finding a belter 
optimization starting point for converging to a lower 
weight design is introduced. Finally, a geometric 
nonlinear analysis is performed to confirm that the final 
design, obtained from linear optimization analyses, 
meets all the design constraints. 

Design Requirements and Loads 

All the tank structures arc assumed to be symmetric 
sandwich constructions. One of the most important 
design requirements for the LOX or LH2 cryogenic 
tanks is that the minimum tank wall thickness can meet 
permeability constraints. For the sandwich lank wall, 
the minimum inner facesheet thickness is 0.06 in. and 
the strain allowable of the facesheet is 6,000 
microstrain (6.0x10 ^ in./in.). However, further 

analyses and experimental tests are required to validate 
this design requirement. The thickness of the outer 
facesheet is assumed to he equal to the inner facesheet 
to avoid coupling between bending and extension. 

The design requirements and loading conditions used in 
this tank sizing and analysis study are based on data 
reported in the literature [4, 6]. The present study uses 
the most severe loading condition, which corresponds 
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to a 3-G axial acceleration. The characteristics of this 
loading condition are listed in Table I. 


Table 1. Maximum 3-G Axial Acceleration Load 
Condition. 


Loading Description 

Value 

Velocity 

4, 140 ft/sec. 

Angle-of-Attack 

-0.044 degrees 

Axial Thrust Force 

4.36x1 OMbs. 

Axial Drag Force 

-1.85x10 s lbs. 

Vehicle Acceleration 

3.0 g’s (Axial) 
0.063 g’s (Normal) 

Tank Ullage Pressure 

25.0 psig (LOX) 

30.0 psig (LH2) 

Fluid Weight 
LOX (46.2% full) 

0.91x10* lbs. 

LH2 (50.0% full) 

0.164x10’ lbs. 


Material Properties 


Material properties representing a quasi-isotropic lay- 
up of toughened graphite epoxy (IM7/977-2) tape 
material are used for sizing the tank structures. Material 
properties and minimum gage requirements are 
presented in Table 2. 

Table 2. Material properties and minimum gages. 


Property Description Value 


Quasi-Isotropic IM7/977-2 Laminate 
Young's Modulus (Msi) 8.804 

Shear Modulus (Msi) 3.354 

Poisson’s Ratio 0.312 

Density (Ib./in. 1 ) 0.057 

Tank wall strain allowable 0.006 

(lension/compression, in. /in.) 

Minimum Gage* - Tanks ( in.) 0.06 

Minimum Gage* - Intertank (in.) 0.04 


Minimum Gage* - Internal structures (in.) 0.02 

Honeycomb Core (Phenolic. 3.0 Ib./ft/) 

Shear Modulus (ribbon direction. Ksi) 20.0 

Shear Modulus (transverse ribbon direction, 9.0 

Ksi) 

* Top or bottom sandwich facesheet thickness 

Tank Structure Arrangements and Finite Element 
Models 

The finite element model of the conformal tank with its 
interna! structure is shown in Figure 2. Due to 
symmetry, only half of the lank structure is modeled. 
The ‘exploded view’ of the tank shown in Figure 3 
includes the intertank, the LOX tank, the LH2 tank and 


other internal structures. The intertank is a structure that 
connects the liquid oxygen (LOX) and liquid hydrogen 
(LH2) cryogenic tanks. Internal structures maintain the 
tank shape and prevent large out-of-plane deformations. 
The three translational nodal displacements of the tank 
structure are fixed along the curve ABC near the rear 
end of its LH2 tank, and symmetry boundary conditions 
are applied on the symmetry plane (see Figure 2). The 
aforementioned “Maximum 3-G Axial Acceleration 
Load Condition” is applied to the tank structures. This 
loading condition consists of a transverse load ( in- 
direction load) due to the angle-of-attack. axial load (x- 
direetion load due to axial thrust and drag forces, tank 
ullage pressure, and tank hydrostatic pressure due to 
fluid under vehicle acceleration. The transverse and 
axial loads are applied to the tank outside surface as 
nodal forces. The ullage and hydrostatic pressures are 
applied as distributed pressures on the internal surfaces 
of the LOX and LH2 tanks. 

All lank structures, including the walls of the LH2 and 
LOX tanks, intertanks, and tension ties are honeycomb 
sandwich constructions. Strength of materials formulae 
and the minimum gage requirements are used to set the 
upper and lower bounds of each design variable in the 
optimization analyses. 

Cryogenic Tank Structure Optimization Analysis 

The optimization problem is formulated as 

Minimize weight min VV'(x) 

Subject to constraints 

LOX and LH2 Tank Walls: 
strain 

-6,000 < Principal f(x) < 6.000 microstrain 
total facesheet thickness 
0.12 in.<TF(x) <0.2 5 in. 
sandwich core thickness 
0.25 in. < TC(\) <3.0 in. 

Intertank: 

strain 

-6,000 < Principal £(x) < 6.000 microstrain 
total facesheet thickness 
0. 1 2 in.< TFCx) <0.2 5 in. 
sandwich core thickness 

0.25 in. < TCfx) < 3.0 in. (Cylindrical regions) 

0.6 in. < TC(x) < 3.0 in A Flat regions) 

Internal Structures: 
strain 

-10.000 < Principal f(x) < 1 0.000 microstrain 
total facesheet thickness 
0.04 in<TFix) <0.2 5 in. 
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sandwich core thickness 
0.25 //?.< TC{\) < 3.0 in. 

where X represents the design variables. 

The conformal tank structures are divided into 33 
regions as shown in Figure 4. In each region, the 
composite facesheet and core thicknesses are treated as 
continuous variable. Since each region has two design 
variables, the total number of design variables for the 
conformal tank is 66. The lower bound of the core 
thickness of the intertank Hat regions is set to 0.6 in. to 
ensure sufficient bending rigidity to prevent large out- 
of-plane deformation induced by the compression 
loads. 

The structural optimization of the tank is performed 
using MSC/NASTRANs optimization capability 
(solution sequence number SOL 200) |5|. To perform 
optimization analysis, a design model needs to be 
developed first whose essential features are the 
definition of design variables, the functional relations 
between these design variables and the structural 
sectional properties, the definition of the objective 
function, and the constraint bounds. Detailed 
MSC/NASTRAN input entries for setting a design 
model for honeycomb sandwich structures arc 
presented in the Appendix. To lest the completeness of 
the input entries and check the optimization analysis 
results, optimization analyses of a uniform cross- 
section beam and a three-section beam under a uniform 
distributed load are performed. These results are also 
presented in the Appendix. 

MSC/NASTRAN incorporates three algorithms to 
perform structural optimization, including the Modified 
Method of Feasible Directions (MMFD) [7], Sequential 
Linear Programming (SLP), and Sequential Quadratic 
Programming (SQP) [2,3 j. The Modified Method of 
Feasible Directions is used in most of the optimizations 
perlormed in this study. However, all three algorithms 
are gradient-based optimization procedures and can be 
trapped at a local minimum. Nongradient based 
algorithms such as Simulated Annealing |8] or Genetic 
Algorithms [9], that can prevent the optimizer from 
trapping at a local minimum, are not available in 
MSC/NASTRAN. 

Technical Difficulties with Conformal Tank 
Optimization Analyses 

Two optimization analyses are performed first, one 
starting from the upperbound of ail design variables and 
one starting from the lowerbound of all design 
variables. Weight and maximum strain constraint 
violation vs. cycle number of these two optimization 
analyses are shown in Figures 5 and 6. The maximum 


strain constraint violation is defined as the maximum 
percentage of strain higher than the allowable strain. 
These figures show that optimization analyses can be 
used to reduce the weight of an overly conservative 
design or to modify an infeasible design to become a 
feasible design. The weight of an over-designed 
conformal tank is continuously reduced bv the 
optimizer as shown in Figure 5. Furthermore. Figure 5 
shows that there is no constraint violation at any part of 
the optimization cycle. Figure 6 shows that when 
starting from an infeasible design formed by the 
lowerbound of each design variable, weight is added 
during the initial cycles to create a feasible design. 
After the design becomes feasible, the optimizer then 
lowers the weight of the design. 

The optimal designs obtained by using the upperbond 
and the lowerbound starting points are designated as 
Design A and Design B, respectively. The weight of 
Design A is and 53.983 lbs. and the weight of Design B 
is 36.283 lbs. as shown in Figures 5 and 6, respectively. 
The weight of Design A is about 5(KT heavier than the 
weight of Design B. A study to determine whether there 
are barriers that can prevent the optimizer moving from 
Design A to Design B is conducted. 

In this study, a design line connecting the two optimal 
designs. Design A and Design B. is used. Nine 
intermediate points between the two end design points 
along the design line are created. If the vector of design 
variables ol Design A is jf, and the vector of design 
variables of Design B is . the design points along the 
design line are 

A - = ,-x Aj + ( I - r ) x ,T 2 . r = 0.0,0. 1. 0.2, 0.3 1.0 ( 1 ) 

Linear static finite element analyses are performed to 
compute weight and evaluate constraint violation for 
each design point. It is found that weight is a linear 
junction ol r. However, some of the design points along 
the design line are infeasible designs, which violate 
strain constraints at /-0.7, 0.8 and 0.9 as shown in 
Figure 7. This indicates that the design space is 
noneonvex. Although the maximum strain constrain 
violation is less than 1 .5 r /c. this noneonvex design space 
can create local minima (barriers) to stop optimization 
analyses. For optimization solution convergence, the 
constraint violation tolerance limit is set to ().59f- as 
shown in Figure 7. 

All strain constraint violations occur in an element that 
is located in the LOX tank rear dome (Design Region 5 
in Figure 4) at the middle of the interface line among 
Design Regions 5. 6. and 10 (see Figures 4 and 8). The 
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strain change for this element as a function of r-value is 
plotted in Figure 8. 

Optimization analyses using each design point along 
the design line as a starting point are also performed. 
Two optimization methods, the MMFD and the SLP. 
are applied. The optimal weights are plotted as a 
function of r-value as shown in Figure 9. Note that each 
method predicts a local minimum for each r-value. The 
local minimum weights for each r-value predicted by 
both methods agree very well with each other. A 
maximum difference of about five percent occurs at 
r=0.6. For this tank design there appears to be many 
local minima. A simple procedure that can create a 
better starting point for the optimizer is presented in the 
following section. 

An Improved Optimization Procedure for Weight 
Reduction 

In this procedure, a new starting point is assembled by 
selecting the lowest design value for each design 
variable in any of the aforementioned eleven optimal 
designs obtained by the MMFD. An additional 
optimization analysis is then performed using this new 
vector of design variables as a starting point. The new 
optimal design (Design C) has weight of 35,533 lbs. 
which is 750 lbs. lower than the previous lowest weight 
design (Design B). No further weight reductions of 
significance are found by searching in the vicinity of 
Design C. Furthermore, a few additional optimization 
analyses using various selected starting points that are 
not on the design line are conducted. Optimal weights 
from these analyses are all higher than the weight of 
Design C. 

The design values and strain distributions of each 
region of the final design are then examined. It is found 
that high local strains dictate the design of Design 
Region 21 in Figure 4. Thus, it is expected that splitting 
the high strain elements out to form an additional 
design region may achieve further weight reduction. 
The high stress areas at two corners of Design Region 
21 are isolated to form a new' design region. Thus, the 
Design Region 21 is split into Design Regions 21a and 
21b as shown in Figure 10. An optimization analysis is 
performed with this modified model and a final design 
with a weight of 34,492 lbs is obtained. The design 
variable values and analysis results of this configuration 
are reported in the following sections. 

Thicknesses and Weights of Tank Structures 

The sandwich facesheet thickness, core thicknesses, and 
weight of each design region are listed in Table 3. Most 
of the sandwich tank wall facesheels are at the 


minimum gage thickness. The distributions of total 
facesheet thickness and core depth of the final 
sandwiched conformal tank walls are plotted in Figures 
1 1 and 12, respectively. 

Table 3. Design variable values and weight. 

Design Total Facesheet 


Core 

Weight 

Region Thickness (in.) 

Thickness (in.) 

(lbs.) 

1 

0.12 


0.25 

935 

2 

0.12 


0.25 

681 

3 

0.12 


0.28 

687 

4 

0.17 


0.66 

1283 

5 

0.19 


0.45 

1489 

6 

0.23 


0.34 

844 

7 

0.12 


0.25 

33 

8 

0.12 


0.26 

98 

9 

0.15 


0.76 

226 

10 

0.04 


0.28 

579 

1 1 

0.04 


0.28 

501 

12 

0.04 


0.28 

282 

13 

0.04 


0.27 

100 

14 

0.04 


0.27 

1 (X) 

15 

0.12 


1 .30 

1 189 

16 

0.12 


0.30 

2169 

17 

0.12 


0.28 

945 

18 

0.12 


0.34 

1047 

19 

0.12 


0.82 

1 164 

20 

0.15 


1.63 

1940 

21a 

0.15 


0.29 

1060 

21b 

0.25 


1.54 

83 

22 

0.12 


0.31 

1099 

23 

0.15 


0.61 

902 

24 

0.12 


0.44 

867 

25 

0.14 


2.83 

1695 

26 

0.25 


3.00 

4288 

27 

0.04 


1 .70 

1213 

28 

0.09 


0.51 

1561 

29 

0.09 


0.79 

351 

30 

0.06 


0.27 

2424 

31 

0.04 


1.61 

876 

32 

0.04 


0.26 

884 

33 

0.04 


0.26 

897 


Total Weight: 34,492 lbs. 


The weight breakdown of the five major components of 
the final design is shown in Figure 13. Note that a 
significant portion of the weight, about 29% of the 
conformal tank weight, is due to LOX and LH2 tank 
internal structures. Therefore, potential further weight 
reduction may be achieved by changing the design of 
the internal structures. 
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Linear and Nonlinear Structural Analysis Results 


Concluding Remarks 


The MSC/NASTRAN design optimization procedure 
computes design sensitivities based on linear static 
finite element analysis results. Geometrically nonlinear 
eflects are not considered. Thus, geometrically 
nonlinear analyses need to be performed to confirm that 
the final tank structure can meet the design 
requirements. If the nonlinear analysis results indicate 
some regions violate the design requirements, these 
regions need to be resized. 

Nonlinear analyses of the final design of the conformal 
lank are performed. An initial nonlinear analysis reveals 
that a deformation larger than that expected occurs in 
the Hat portions of the intertank; thus the core thickness 
of the inlertank is increased to 2.5 inches. A linear 
analysis and a nonlinear analysis are performed for 
comparison. Plots of major principal strains of the 
outside LOX and LH2 tank walls obtained from both 
analyses are shown in Figures 14 and 15. These results 
show that the linal LOX and LF12 tank- wall designs for 
the conformal lank can meet the design strain allowable 
of 6,000 microstrain. The nonlinear membrane stress 
stiffening effect lowers the strains in the flat regions of 
the conformal LOX tank and the rear cylindrical region 
of the conformal LH2 tank. 

Further Design Improvement Areas 

This study used structural optimization to conduct 
sizing of lank structures in the early design stage. 
Similar optimization procedures can be used in the 
detailed design stage of the tank structures. The 
following three design approaches are recommended 
for further weight reduction. 

1. Local high stress concentration areas within a 
particular design region wall dictate the sandwich 
facesheet and core thicknesses for the entire design 
region. If a design region is separated into a low 
stress design region and a high stress design region, 
significant weight reduction from the low stress 
region is possible. Note that this approach is used 
to reduce the weight of Design Region 2 1 . 

2. More cutouts can be added to the LOX and LH2 
internal structures. In this study only the Design 
Regions 12 and 30 have cutouts. 

3. The nonlinear membrane stiffening effect is 
expected to reduce the bending deformation in 
some tank regions. Using nonlinear analyses 
instead of linear analyses in the optimization 
process may achieve a lower weight design. 


The structural optimization analysis capabilities in 
MSC/NASTRAN were used to size a conformal 
cryogenic tank’s honeycomb sandwich structures. This 
tank is divided into 33 design regions. The facesheet 
thickness and the honeycomb core depth of each design 
region are the design variables. The MSC/NASTRAN 
input entries for setting the design model to optimizing 
these design variables for a minimal weight tank design 
were developed. 

This study found that the feasible design space for this 
problem is nonconvex and many local minima are 
exhibited. Because MSC/NASTRAN does not 
guarantee to find the global minimum, a simple 
procedure, which has potential to create a belter starting 
point for the optimizer, was introduced. This procedure 
was used to obtain a final improved lightweight design. 
Geometrical nonlinear analysis results confirmed that 
this final LOX and LH2 lank-wall designs meet the 
strain constraints. Finally, areas that could have 
potential for further weight reductions were discussed. 
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Appendix 

Design Model for Sandwich Structure Optimization 

MSC/NASTRAN SOL200 input entries to build a 
design model for optimizing sandwich structures are 
developed. A sandwich beam shown in Figure A-l is 
analyzed to test the completeness of these input entries 
and verify the results. The beam is 30 inches long and is 
subjected to a uniformly distributed load of 100 Ibs/in. 
The beam is divided into three different sections. Two 
different optimization analyses are performed: the first 
consists of a uniform beam which has the same 
facesheet and core thicknesses for all three sections and 
the second consists of a three-section beam in which 
each section has different facesheets and core 
thicknesses. The material properties listed in Table 1 
are used. 

The optimization statements are: 

Minimize weight mm Vf(.v) 

Subject to constraints: 
strain 

— 1 OCX) < Axial £{x) < 1,000 micros train 
total facesheet 
QA2in.< TF(x) <0.24 in. 
sandwich core 
0.1 in.<TC(x)< 4.0 i/i. 
where X are design variables. 

The key inputs include the design variables 
(DESVARs), design equations (DEQATNs), design 
variable to analysis property relations (DVPREL1), 
design variable to analysis property relations using 
user-supplied DEQATNs (DVPREL2), design 
sensitivity response quantities (DRESPl), design 
constraints (DSCONs) and design optimization 
parameters (DOPTPRM). 

The DEQATNs inputs related to MSC/NASTRAN's 
PSHELL data entry [ I Oj are 


DEQATN 101 
DEQATN 201 
DEQATN 301 
DEQATN 401 
DEQATN 501 


NSM(TFJC) = TC*RHO 
Zl(TFJC) = -0.5*(TC+TF) 
Z2(TFJC) =0.5*(TC+TF) 
ST(TFJC) = TC/TF 
Bl(TFJC) = 


3. 0 *(TC/TF+0. 5) * *2+0. 25 
where TF is the total facesheet thickness (top 


and bottom facesheets). 


TC is the core thickness, 


RHO is the core density. 

NSM is the honeycomb core areal density, 
Z 1 and Z2 are the fiber distances from the 


reference surface, 

ST is the transverse shear thickness 
ratio (transverse shear uniformly 
carried by the core only), 
and HI is the bending moment initial ratio: 


Bl = 


12x7 
TF x 


(A- 1 ) 


where the bending rigidity (/) of the beam can be 
expressed as 


/ 


(TF)' 

48 


TF > 
TF(TC + — )' 

i 


(A-2) 


The DVPREL1 and DVPREL2 inputs are created based 
on the MSC/NASTRAN's shell element property 
(PSHELL) format using DESVARs and DEQATNs. 

The DRESPl entries define a set of structural responses 
that are used in the design either as constraints or as an 
objective. The DSCONS entries define design 
constraints in design sensitivity analysis. The 
DOPTPRM entries are used to override default values 
of parameters used in design optimization. 

Beam Results 


The optimized uniform cross-section beam has a 
facesheet thickness of 0.159 in. and a core thickness of 
3.047 inches. In the optimized 3-section beam. Sections 
1 and 3 have the same facesheet thickness of 0.146 in. 
and same core thickness of 2.742 in; and Section 2 has 
a facesheet thickness of 0.153 in. and a core thickness 
of 3.173 in. The weight for the uniform cross-section 
beam and the three-section beam are 0.408 and 0.383 
lbs, respectively. 

The top-surface axial strain plots of the uniform beam 
and the 3-section beam are shown in Figure A-2. The 
strain results of the uniform beam agree with the closed 
form solutions. The 3-section beam has a much longer 
section subjected to a high strain near the maximum 
strain of 1.000 microstrain as expected. 
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Tension Ties 

Figure 3. Substructures of conformal tank configuration 
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Figure 5. Optimization analysis starling lroni the upperbound of all design variables. 
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Figure 6. Optimization analysis starting from the lowerbound of all design variables. 
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Figure 7. Infeasible design at r=0.7, 0.8, and 0.9. 



Figure 8. Strain changes as a function of r-value at the constraint violation location. 
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Figure 1 3. Weight distributions of major tank substructures. 
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Figure 14. Strain plot of conformal tank from linear analysis. 
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Figure 15. Strain plot of conformal tank from nonlinear analysis. 
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Figure A-l. Beam under a uniform distributed load. 



Axial Strain. 
xKT 6 in./in 

0 

-950 | 

-1000 * 


Figure A-2. Axial strain plot for the beam top surface. 
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